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1. Introduction 

Transport coefficients characterize fluctuations and relaxation on long length and time scales 
in systems slightly away from thermal equilibrium. Consider a conserved current, dfU + V • j = 0, 
in a system with a small nonuniformity in the density n{t,x). To first approximation, a current will 
form to wash out the imbalance: j = —DWn + . . ., where the dots indicate higher order terms in 
gradients. Combining the conservation law with this constitutive equation immediately yields the 
diffusion equation dtn = DW^n. The proportionality factor D is the diffusion constant and is the 
only quantity determined by the underlying microscopic theory: one may view it as a low-energy 
constant. Similarly, consider the conserved energy-momentum tensor T^^ . While in equilibrium 
jij — §ijp^ vvith P the pressure, a perturbation characterized by a nonuniform flow field u{t,x) will 
change this to (in the local rest frame, where r"' = 0) 

T'j = d'jp - T] (^d'uJ + dJ'u' - - i;d'jd,u' + .... (1.1) 

Again the dots indicate higher order terms in gradients. The coefficient of the traceless combina- 
tion is the shear viscosity Tj and is the bulk viscosity. Combining energy-momentum conserva- 
tion with the constitutive equation above yields the hydrodynamic equations, predicting e.g. sound 
waves. In the context of thermal QCD one may therefore view hydrodynamics as the low-energy 
effective theory describing real-time dynamics at sufficiently large length and time scales. The 
form of the hydrodynamic equations is fixed by combining exact conservation laws and constitu- 
tive equations, which are obtained in a gradient expansion. In the latter, a number of low-energy 
constants appear, determined by QCD: shear viscosity rj, bulk viscosity i^, electrical conductivity 
a, diffusion constants D, etc. In this talk I review the progress in determining these coefficients 
from first principles. 

The recent interest in transport in QCD and related theories is mainly due to the relativistic 
heavy ion program at RHIC The remarkable effectiveness of ideal hydrodynamics in describing 
heavy ion phenomenology ^ suggests that transport coefficients are very small (when appro- 
priately normalized). This in turns implies that thermalization times are short and interactions are 
strong, suggesting that in the temperature range 1 ^ T/Jc < 2 the quark-gluon plasma is not a 
weakly coupled system of quarks and gluons, but instead strongly interacting (sQGP). A second 
reason for interest in transport is the fertile applicability of gauge/gravity duality (or AdS/CFT 
correspondence) to study strongly coupled thermal gauge theories in the hydrodynamic regime 
[Q]. This has led to many (semi-)analytical results for those theories that admit a gravity dual and 
provides an important stimulus for QCD, where a gravity dual is not available. The best-known 
example concerns the ratio of the shear viscosity and the entropy density s, 

n I 
s An 

which is obtained in all thermal gauge theories in the (strongly coupled) regime described by a dual 
gravity theory. This ratio is much smaller than in weakly coupled QCD, where 

TJ 1 

lim — ~ > oo. (1.3) 

g^o s g^lnl/g 

It is an open question what this ratio is in QCD just above the deconfinement transition. Below I 
describe recent progress in lattice QCD that will bring us closer to answering that question. 
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2. Transport coefficients at weak and strong coupling 

According to the Kubo relations, transport coefficients are given by the slope of current-current 
spectral functions, computed in thermal equilibrium, at vanishing energy. For example, the electri- 
cal conductivity a (or charge diffusion constant D times susceptibility S) is determined by 

a = DS=lim ^ ^ \ (2.1) 
cB^o 2(0 

where the current-current spectral function is given by 

p^'ico) = J d'xe''''{[fit,x),fm)e^- (2.2) 

For a single fermionic charge carrier with charge e the electromagnetic current reads = exj/y^Y- 
Similarly, the shear and bulk viscosities are given by 

77 = lim C = :^ lim ^-^A^, (2.3) 

' ffl^o 2(0 ^ 9 £0^0 2(0 

in terms of the energy-momentum tensor spectral function 

pMv,pa(^) ^ J d\e""{[T^\t,x)JP%0)]),^. (2.4) 

2.1 Weak coupling 

At weak coupling transport coefficients can be computed using either kinetic theory or by 
summing sets of Feynman diagrams. In this limit, transport coefficients are manifestly proportional 
to the mean free path or the inverse collisional width l/F. They are therefore very large and 



inversely proportional to the coupling constants in the theory, see Eq. (1.3). In ultrarelativistic QCD 
the shear viscosity, electrical conductivity and diffusion constants have been computed to leading- 
logarithm order ^ and subsequently to full leading order ^ in the gauge coupling. The relevant 
physics at leading-log are those 2 •s-^ 2 scattering processes that are logarithmically sensitive to 
infrared screening effects. The extension to full leading order requires, besides the other 2 2 
scattering processes, also the inclusion of specific particle number changing processes [^. The 
bulk viscosity is more complicated since particle number changing processes need to be included 
from the start; the first calculation in QCD can be found in Ref. [||]. It is shown that a parametrically 
correct estimate is given by 

C ~ 1577 (1/3 -V?)', (2.5) 

where for light quarks the speed of sound is determined by — 1/3 ~ j8(g^) ~ g^. The equiv- 
alence between kinetic theory and diagram summation in thermal gauge theories has been demon- 



strated in a number of papers [Q, |10[ [11|, [12|, |13|]. A one-loop calculation is never sufficient; instead 
an infinite set of ladder diagrams has to be taken into account. 

A complete leading-order calculation is also possible in the large Nf limit, where Nf indicates 
the number of flavours. This has been done using kinetic theory for massless quarks |[T^] and 



extended to massive quarks using diagrams [15]. Only Coulomb scattering processes contribute. 
Also in this limit one is effectively in the weakly coupled regime, since r\ / s ^ Nf ^ °°. 
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Figure 1: Typical skeleton diagram contributing to the shear viscosity in the 0{N) model at large A^. 



In all cases above only the leading-order result is known. This prevents any reasonable extrap- 
olation to stronger coupling, unlike in the case of e.g. the pressure where the terms up to g^ln 1/g 
are known [|T^. Recently the first subleading correction to a transport coefficient in a relativistic 
field theory has been computed by Moore [17], namely to the shear viscosity in scalar theory. 
The shear and bulk viscosity for single-component scalar theories have been computed a long time 
ago to leading order at weak coupling using ladder diagram summation [ITSl]. In the 0{N) model 



the shear viscosity has also been found to leading order in the \/N expansion [|l9|], using similar 



techniques as in large Nf gauge theories. It is shown in Ref. [ ]17[ ] that the first subleading correction 
to the shear viscosity at weak coupling is sensitive to soft physics and can be extracted from the 
same set of diagrams that contribute at leading order by an expansion in terms of the thermal mass 
mth/r ~ The result in the 0{N) model, with & X / {A\N){<^a<^af interaction, reads [|l7]| 



3033.54 + 1548.3^ 



3033.54+1548.3 



+ 



(2.6) 



where the second line is valid in the combined large N and weak coupling limit. As mentioned 
above, the complete large A'^ result, without employing the weak-coupling expansion, is also known 



[19]. This includes a resummation of the thermal mass to all orders, as well as the inclusion of 



other diagrams suppressed in the weak coupling limit, see Fig. Therefore, an assessment of the 
perturbative expansion can now be made. This is shown in Fig. ^ where the perturbative result 




Am 



Figure 2: Shear viscosity Tj, normalized with the result at vanishing coupling, in the 0{N) model as a 
function of \^X: comparison between the perturbative [|l7|| and the exact [|l9[ result in the large limit. 
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T] /t]o = 1 + is compai-ed with the complete large N result. Here T]o is the result at vanishing 
coupling, which is used for normalization. From the first two terms, it appears that the convergence 
of the weak coupling expansion is not impressive. A possible approach to improve the convergence 



[ ]17[ ] would be to include the thermal mass without re-expanding, as has been tried for the pressure 
[20]. How to include all ^(A) effects is currently not known. 

Considerable effort has been spent on nonrelativistic dynamics and diffusion of heavy quarks 
(with Mq » T). I will not discuss that here, but refer to Ref. [21 1 where definitions and a leading- 



order calculation of the momentum diffusion coefficient can be found, and to Ref. [ |22[ ] for the first 
next-to-leading order result at weak coupling. 

2.2 Strong coupling 

In order to complement weak coupling results as discussed above, it is desirable to have an 
analytic method tailored for strongly coupled field theories. For a certain class of thermal gauge 



theories such a method is available, via the gauge/gravity or AdS/CFT correspondence [|23|]. It uses 
the duality between a field theory at finite temperature and black holes in AdS space described in 
string theory. The best-known case relates = 4 supersymmetric Yang-Mills theory to type IIB 
string theory on AdSs x , but there are many other field theories for which a gravity dual can be 
found (but not QCD). The parameters in field theory, the number of colours and the coupling 
constant g^, are related to the parameters of string theory in AdS space. It turns out that the duality 
is most powerful when is small, since then loops in the string theory are suppressed, but with the 
't Hooft coupling A = g^Nc large, since then stringy effects are suppressed and the string theory 
reduces to a supergravity theory. Strongly coupled gauge theories, with large 't Hooft coupling, 
are therefore the natural ai^ea of application. Using the gauge/gravity correspondence, it has been 
shown that the field theories for which the duality holds behave hydrodynamically at large length 
and time scales, which supports the framework of (nearly ideal) hydrodynamics to understand the 
dynamics in thermal gauge theories at strong coupling. 

The most famous example concerns the shear viscosity in ^ = 4 SYM. In the Unlit that both 



Nc and A go to infinity, it is equal to [ 24 j 



^N^T^. (2.7) 



In the same limit the entropy density is [^] 



J ^ 3 
—N^T^ = -s 
2 ' 4 



such that the ratio is Tj/s = 1 /4n. This result for the ratio is universal and is achieved in all thermal 
gauge theories that can be described by a gravity dual [26, |^]. The bulk viscosity vanishes due to 
conformal invariance. 

The shear viscosity in ^ = 4 SYM at strong coupling is another rare example of a transport 



coefficient where the first subleading correction is known. One finds [ )28p 

75C(3: 
4AV2 
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The first correction to the ratio then reads 



^ 1 / 135C(3) 

- = 1 + . + • • • • (2-9) 



471 V 8A3/2 

It has been conjectured that 1/471 is a lower bound for a wide class of systems, including e.g. 



the quark-gluon plasma and trapped atomic gases [ |29| , [Zq ]. For a recent critical assessment, see 



Ref. [30] 



Other transport coefficients that have been computed include the /?-charge diffusion coefficient 
and the bulk viscosity. These are not universal and depend on the theory under consideration. In 
=yK = 4 SYM the /?-charge diffusion constant, susceptibility and conductivity are given by [ ^l] ] 

1 A^2j2 

D = - , S = ^— , <^ = -^- (2.10) 

2-nT' 8 ' 1671 ^ ' 

The bulk viscosity in a number of theories is discussed in Ref. [p^]. For strongly coupled systems 



it is proposed that it satisfies (cf. Eq. ( [2. 51 ) and the absent square) 

C>2t](1/3-v2). (2.11) 
Also heavy quark dynamics has been studied extensively in strongly coupled ^ = 4 SYM, see e.g. 



Ref. [ |33| ] for a clear discussion. 

Due to the interest in .yV = 4 SYM as a testbed for thermal field dynamics, various quantities 
have also been studied at weak coupling. The ratio r]/s has been computed in this theory to full 
leading order in the small A limit (the result is independent of Nc) and compared with the same 
quantity in QCD [Q]. Heavy quark diffusion has been studied to leading order at weak coupling 
in Ref. 



3. Spectral functions 

So far, in discussing transport coefficients, the focus has been on the zero-energy limit of 
spectral functions. In this section I give a (restricted) review of spectral functions at arbitrary (o, 
motivated mainly by hydrodynamic structure at weak and strong coupling. This will turn out to be 
useful for lattice QCD studies of spectral functions, to be discussed next. 

In the weak coupling limit, current-current spectral functions have a characteristic energy de- 
pendence [^. At large energy they increase as ft)", where the power n is determined by the (mass) 
dimension. For example, the energy-momentum tensor spectral function (2.4) increases as (O^ and 
the EM current spectral function (2.2) as g)^ (unless there is a cancelation between components). 



At small energies, there is a transport peak. In free field theory, this peak manifests itself as a 
singular term, 

^~27r5M, (3.1) 

CO 

reflecting that in a free theory the mean free path is infinite and transport coefficients diverge. 
Interactions regulate this singular behaviour and after the resummation of the collisional width F, 
the transport peak is modified to 

m ^ (32) 
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Figure 3: Shear viscosity: p{(o)/(oT^ vs. co/T in QCD (sketched) for two values of the colhsional width 
r/T = 0.1,0.5 (left) and in ^ = 4 SYM at strong coupling (right). 




Figure 4: Charge diffusion: p{co)/coT vs. co/T in QCD (sketched) for two values of the collisional width 
r/T = 0.1,0.5 (left) and in ^ = 4 SYM at strong coupling (right). 



For simplicity, I parametrized the effect of interactions with a single constant width F; in reality 
the width depends on the momentum of the (quasi-)particles contributing to transport. At weak 
coupling the collisional width F g^T: therefore the transport peak is narrow {r^ g^) and high 

The (o dependence of spectral functions of conserved charges, such as = j^T'^^{t,iL) or 
Q = J^f{t,x), is completely fixed by the conservation laws. For example, the spectral function for 
the total charge density reads 

^-^=Z27i5{(0), (3.3) 

(0 

where S is again the susceptibility. This ensures that unequal-time correlation functions, such 
{Q{t)Q{0)), are in fact constant and equal to VTE, where V is the spatial volume and it is assumed 
that {Q) =0. Here I only discuss spectral functions at zero spatial momentum; for nonzero spatial 
momentum see e.g. Ref. [37|. 

Combining the rising high-energy part and the transport peak at small energies yields spectral 
functions as shown in Fig. ^ (left) for the energy-momentum tensor (in the channel relevant for the 
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shear viscosity) and Fig. ^ (left) for the EM current (charge diffusion). Here the width is taken 
as a free parameter and artificially increased to larger values (T /T < 1), while keeping the rest of 
the spectral function unchanged. As a result the transport peak, which is narrow at weak coupling 
(F/r <^ 1), becomes less singular and gets smeared out. 

Interestingly enough, this behaviour can be compared with spectral functions computed in 



jy = 4 SYM in the limit of large Nc and large 't Hooft coupling A, using gauge/gravity duality [38 



. The energy-momentum tensor spectral function has to be computed numerically by solving an 



ordinary differential equation, but for the /?-current the analytical result is known and reads [40] 



_ N'^ CJ^sinh(o)/2r) 
P^^>~ 16a;cosh(w/2r)-cos(w/2r)' ^ ' 

These spectral functions are shown in Figs. ^ and || (right). Rather than dividing by powers of A^^c. I 
have simply put Nc = 3. It is clear that at strong coupling the transport peak is no longer separated 
from the high-energy contribution and the spectral functions go smoothly to ft) = 0. The intercept 
is of course proportional to the shear viscosity and the conductivity respectively. 

Comparing the spectral functions on the left with those on the right makes it interesting to 
speculate what happens in QCD in the strong coupling regime above the deconfinement transition. 

4. Transport from lattice QCD 

We have seen that weak-coupling methods are probably not applicable in the interesting tem- 
perature regime of QCD probed by current heavy ion collisions. Furthermore, these calculations 
are so involved that in most cases only the leading-order result is currently known; it is an open 
question how to determine subleading corrections. On the other hand, the strong-coupling results 
discussed above are obtained in theories that are not QCD. So the important question is what can 
be said about transport in QCD when 1 ^ r/Jc < 3, using nonperturbative lattice simulations. 

The euclidean correlator calculated in numerical simulations is related to the corresponding 
spectral function via a dispersion relation, 

G,(,«*,)=r^5^-^^. (41) 

J-oo IK (O — 1(0,, 

where ft)„ = 2nnT {n G Z) are the bosonic Matsubara frequencies. In eucUdean time, with < T < 
1/r, this relation becomes 

f°° d(0 

Ge{t)= / —K{(0,T)p{(0), (4.2) 
Jo 271 

with the kernel 

= = [1 + + (4.3) 

smh(ftj/2rj 

where nB{(o) = Xjie^l^ — 1) is the Bose distribution. The first expression for the kernel shows 
the characteristic euclidean time dependence, while the second expression emphasizes that the 
correlator and its spectral function are essentially related via a Laplace transform, made periodic to 
satisfy the Kubo-Martin-Schwinger (periodicity) condition. 
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If the euclidean correlator is known analytically, the spectral function can simply be obtained 
by analytic continuation, 

p{(0) = 2\TaGE{i(0n (0 + ie). (4.4) 

However, if G£(t) is determined numerically on a lattice with A^^ = 1 /aT points in the euclidean 
time direction (a is the temporal lattice spacing), Eq. ( |4.1D has to be inverted by some other means. 
This is an ill-posed inversion problem, since Ge{t) is known at, say, ^(10) data points, whereas 
P(g)) is needed at ff{lO^) values (after imposing a high-energy cutoff OVax and discretizing the 
resulting finite interval < (O < ODmax)- 

One possible solution is to provide an Ansatz for the spectral function, with a small number of 
free parameters. In this case it is important to be able to judge the applicability of the Ansatz. An 
orthogonal approach is to avoid giving functional forms but only supply a minimal amount of prior 
information, such as positivity {o)p{(o) > 0) and asymptotic behaviour (p(w) ~ (O" for large ft)). 
Methods based on this approach are usually collectively referred to as Bayesian techniques. 

For a weakly coupled quark-gluon plasma, the extraction of transport coefficients is notori- 
ously difficult, since euclidean correlators are remarkably insensitive to the structure of spectral 



functions at energies (O <^ T ||3g, |41|]. However, at stronger coupling the transport peak is much 
broader and the small co limit is no longer singular As discussed above, smooth spectral functions 
are also found in ^ = 4 SYM at strong coupling. This opens up the possibility that transport 
coefficients are accessible in lattice QCD above the deconfinement transition. 

So far the number of papers in the literature that have attempted to extract transport coeffi- 
cients from the lattice in a head-on approach is very small. A first attempt to measure transport 



coefficients can be found in the pioneering paper by Karsch and Wyld, using an Ansatz |42]. This 



method was followed by Nakamura and Sakai for the shear and bulk viscosity [43]. For a critical 



discussion of the Ansatz, see Ref. [|36|]. S. Gupta used Bayesian methods to isolate the transport 
contribution at small energies in the case of the electrical conductivity Q]. In the past six months, 
two significant steps have been made. A standard approach to perform the analytical continuation 
using Bayesian techniques is known as the Maximum Entropy Method (MEM), to be discussed 
below. It was known from previous work that MEM performs poorly at small energies. Aarts, 
AUton, Foley, Hands and Kim have identified and resolved a numerical instability in MEM in the 
limit that ftj — > and applied the new formulation to obtain the electrical conductivity ||45|]. Pre- 
cisely determined correlators are essential to have control over the analytic continuation; while for 
the electrical conductivity this is not a problem, for the shear and bulk viscosity standard mea- 
surement techniques are insufficient. Meyer has applied a two-level algorithm to better determine 



energy-momentum correlators and found a result for the shear viscosity [146|]. All calculations to 
date have been performed in quenched QCD, so one may think of the electrical conductivity in 
pure gauge theory as representing the transport properties of a single electrically charged quark 
diffusing through a gluon plasma. 

In the following I describe in some detail the Maximum Entropy Method and indicate why 



the standard algorithm is unstable in the small CO region |45]. In MEM one reconstructs the most 
probable spectral function by extremizing the probability distribution P[p\GH], i.e. the probability 
to find p, given the correlator G and prior information H. Using an identity for conditional prob- 
abilities, P[pG\H] = P[p\GH]P[G\H] = P[G\pH]P[p\H], P[p\GH] is written as the product of a 
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Figure 5: First four basis functions M,(a)) in the singular value decomposition as a function of ao for 
aftfrnax = 5, No, = 1000, Nr = 24, using the standard (left) and the redefined kernel (right). The inset shows 
a blow-up of the small energy region. 



standard likelihood function, P[G\pH] ~ e^^ {x^ fit), and a prior probability, P[p|//] ~ e"'^ which 
is independent of the data. P[G\H] is a normalization factor. The prior information is encoded in 
the entropy term, 



p{(0) -m{(0) -p(co)log- 



(4.5) 



'm(G)) 

via the default model m{(o). Therefore, the combined function to extremize is P[p\GH] ~ ^-L+as^ 
where a determines the relative weight of the data versus the prior information [47]. 

The most important aspect for my purpose here is the reduction step. Recall that after dis- 
cretization p{(o) is wanted at A^o = ^(10^) values whereas the correlator is only known at ^(10) 
points. I denote with the number of time slices included in the analysis; due to reflection sym- 
metry, N < Nt/2. To make this well-defined, the number of coefficients parametrizing the spectral 
function cannot exceed N, or in other words, p{(o) has to be restricted to an (at most) N dimen- 
sional subspace. In Bryan's algorithm [48] this is achieved via a singular value decomposition 
(SVD) of the kernel K{(On, f,). Viewed as an NaxN matrix, the kernel is written as K = UWV^ , 
where U is an Na x N matrix, with U^U = ^nxN, W is a diagonal N x N matrix, and V is an 
orthogonal N x N matrix. The N dimensional subspace is spanned by the column vectors of U : 
m, (g)„) = Uni- These basis vectors are orthogonal but not complete. An analysis of the extremum 
conditions shows that it is natural to write the spectral function in terms of these basis vectors 
as p{(o) = m(ft))exp^^jC,M,(ft)). This ensures positivity and provides the actual reduction step. 
Extremizing the probability distribution leads to nonlinear equations for the A'^ coefficients c,. 

The first four basis functions in the SVD are shown in Fig. ^ (left), for a typical choice of 
A'^ = A^t/2, A^o, and COmax- A blow-up of the small energy region reveals that the basis functions 
appear to diverge when ft) — > 0, although they are normalized (the smallest energy included here is 
aAco = a(Omax/Nco = 0.005). This apparent divergence is due to the singular behaviour of the kernel 
(p3|): in the limit that ftj ^ one finds that K{(0,t) = 2T/(0 + ^{(o/T). Note that the leading 
singular term is T independent; all T independence resides in the subleading terms. The behaviour 
of the basis functions is therefore indicative of a real problem, which cannot be solved by e.g. 
decreasing Aco. In actual applications of MEM, we (and others) found irregular behaviour at small 



10 



Transport and spectral functions in high-temperature QCD 




_. N =1000, b=Q.l 

ol , \ , \ , \ ^ , 1 

2 4 6 8 10 

co/r 

Figure 6: Current-current spectral function p {(£>) / coT, where / = ifrYy/, as a function of co /T in quenched 
QCD for T/Tc ^ 1 .5 (hot, A^^ = 24) and 2.25 (very hot, A^^ = 16). The intercept at O) = is proportional to 
the electrical conductivity. Results are shown for Na = 1000, 2000 and b= 1.0, 0.5, 0.1 at fixed afiJmax = 5. 



ft), in particular at the smallest nonzero value of ft). Needless to say that this prevents access to 
the transport properties encoded in euclidean correlators. It is worth pointing out that this problem 
only appears at finite temperature, since at zero temperature the kernel reduces to K{co, t) = e^®^ 
and the limit ft) ^ is smooth. It is a manifestation of the fact that the limits (O —^0 and T ^ do 
not commute, which is well-known in thermal field theory. 

Once the problem is identified, it is straightforward to solve it. The 1 /ft) divergence can be 
avoided by defining 

K{o),t) = —K{(o,t), p[co) = —p{(o). (4.6) 

21 (0 

Since K{o), T)p (ftj) = K{o), T)p (ftj) the standard relation with the euclidean correlator holds. How- 
ever, the modified kernel is finite when co ^0: K{0, t) = 1. A SVD of K yields new basis functions 
Ui{co). The first four are shown in Fig. ^ (right). Clearly they take a finite value when ftj — > 0, and 
the (0 = point can be included in the analysis. The redefined spectral function is expanded as 
p(ftj) = rn{(o) expY^LiCiUiio}) and the same MEM routine can be used to find the coefficients c,. 
MEM now reconstructs p ~ p/ftj rather than p. This reshuffling of powers of CO is nontrivial, since 
p and p are not expanded in a complete set: the reduction step restricts p to a different subspace 
with manifestly different properties, in particular in the small CO limit. 

A second (minor) modification is needed to access p{co)/co at zero CO, relevant for transport 
coefficients. For spectral functions of fermion bilinears, such as = yj^Y' the traditional default 
model is m{co) ^ m{co)/co ~ CO, determined by the high-energy behaviour p(ftj) ~ co^. Unfortu- 
nately, this introduces a bias and puts the intercept equal to zero from the start. To avoid this, one 
may use m(ftj) ~ (ft + co), where ft > is a parameter that can be used to assess default model 
dependence at small CO. 

We have applied the modified algorithm to the problem of the electrical conductivity (or charge 
diffusion) in quenched QCD with light staggered fermions and performed simulations on a fine lat- 
tice at j8 = 7.192 of size 64^ x A^^ above the deconfinement transition. The quarks are so light that 
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chiral symmetry restoration is clearly visible when comparing pseudoscalar and scalar correlators 
[p9|]. In Fig. 1^ the spectral function normalized by coT is shown for not too large CO at two tempera- 
tures above the deconfinement transition: T/r^ ~ 1 .5 {N^ = 24) and T jT^^ 2.25 (A^t = 16). In both 
cases results for three values of the default model parameter h are shown. The slight spread of the 
curves gives an indication of the uncertainty in the MEM reconstruction. There is no dependence 
on the discretization along the ft) axis, as can be seen from the results at 'Nx = 24 using A^o = 1000 
and 2000 at fixed ftViax- From the intercept the conductivity is found to be a/r = 0.4 it 0.1 with 
no significant temperature dependence. This result is normalized to a single flavour and should be 
multiplied with the sum of the electric charge squared for light flavours. The error is systematic and 
due to the MEM uncertainty. The statistical error is expected to be smaller. This result is indicative 
of strong interactions: at weak coupling the conductivity behaves as a/r ~ Xj (g'XxvXj g) — > oo, 
whereas at strong coupling the scale is set solely by the temperature. 

A similar conclusion has been drawn by Meyer in the case of the shear viscosity. In Ref. [46] 
an upper bound T]/^' < 1 is obtained in SU(3) gauge theory on lattices at T jTc = 1-24 (j8 = 6.2, 
20^ X 8) and 1.65 (j3 = 6.408, 28^ x 8). 



5. Quarkonium at high temperature 

Another signal of strong interactions in the quark-gluon plasma is the survival of charmonium 
and other heavy quark mesons above r^.. This has been studied on the lattice a few years ago using 



the Maximal Entropy Method [|50|, gjj, |52|]. A more recent extensive study of charmonium and 



bottonium spectral functions can be found in Ref. [p3]. This topic has been discussed last year in 



Hatsuda's plenary talk [|54|], so I briefly mention developments that took place recently and were 
discussed at this Conference. 

Up to last year, all studies were performed in quenched QCD. The TrinLat collaboration has 
carried out dynamical simulations with two flavours on highly anisotropic lattices, which can be 
used to study charmonium at zero and nonzero temperature [^]. A spectral function analysis above 
Tc, using the modified version of MEM described in the previous section, can be found in Ref. [p^]. 
The results suggest that the S-waves {J and T]c) survive up to temperatures close to 27^, while 
the P-waves {Xm and Xcx) melt away below 1.27^. However, there are systematic uncertainties that 
need to be improved in order to make these conclusions more firm, in particular simulations at a 
finer lattice spacing would be desirable. One reason this is necessary is to better understand the 
appearance of artefacts at larger energy introduced by the finite lattice spacing, which have been 
discussed for a variety of lattice fermion formulations (Wilson, staggered, domain wall, overlap) in 
the free field limit [|7|, 1^, p]]. 



The presence of the transport contribution at small « can interfere with spectral features at 



larger O) when not properly disentangled, as emphasized by Umeda [59]. This can be partially 
avoided by subtracting the midpoint value at T = l/2r from the correlator, which has the effect of 
suppressing the contribution at smaU o). A spectral analysis can then be done on G(t) — G( 1 jlT). 

Petreczky and Mocsy have provided a closer look at potential models, traditionally used to 
study quarkonium at zero temperature. Doubts whether potential models can describe quarkonium 



correlators at finite temperature have been expressed in Ref. [pOl]. A lower melting temperature than 
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usually obtained with MEM is found when a potential model based on lattice QCD simulations is 
used, both in the quenched approximation [^] as well as in the theory with 2 + 1 flavours |62]. 

In order to avoid conceptual problems with the extension of zero-temperature potential mod- 
els to nonzero temperature, a real-time static potential, firmly based in thermal field theory, was 
introduced in Ref. [63]. When applied to quarkonium ||6^], the results seem to support the standard 
interpretation of results obtained with MEM. 



6. Summary 

Results from relativistic heavy ion collisions at RHIC have highlighted the importance of 
understanding transport and hydrodynamical behaviour in QCD above the deconfined transition. 
Nonperturbative first-principle calculations of spectral functions, especially at small energies CO < 
T, are badly needed. Since this involves inherently real-time physics, it is a difficult problem for 
lattice QCD, but recently several steps forward have been made. Using multi-level algorithms, 
accurately determined euclidean correlators of the energy-momentum tensor are now available. 
Concerning the analytical continuation to real time, an instability at small energies in the standard 
Maximum Entropy Method, preventing access to hydrodynamical features of spectral functions, 
has been found and resolved. The first results support the idea that the quark-gluon plasma is 
strongly interacting in the temperature range 1 < T /T^- < 2.5. Extension of the work described 
here will hopefully yield a better understanding of the hydrodynamical regime of thermal QCD 
from first principles. 
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